Measurement of the atmospheric neutrino-induced upgoing muon flux using MACRO 
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We present a measurement of the flux of neutrino-induced upgoing muons (<E„ >~ 100 GeV) 
using the MACRO detector. The ratio of the number of observed to expected events integrated over 
all zenith angles is 0.74 ±0.036(stat) ±0.046(systematic) ±0.13(theoretical). The observed zenith 
distribution for —TO < cosf) < —0.1 does not fit well with the no oscillation expectation, giving a 
maximum probability for % 2 or 0.1%. The acceptance of the detector has been extensively studied 
using downgoing muons, independent analyses and Monte-Carlo simulations. The other systematic 
uncertainties cannot be the source of the discrepancies between the data and expectations. 

We have investigated whether the observed number of events and the shape of the zenith dis- 
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tribution can be explained by a neutrino oscillation hypothesis. Fitting either the flux or zenith 
distribution independently yields mixing parameters of sin 2 26 = 1.0 and Am 2 of a few times 10 -3 
eV 2 . However, the observed zenith distribution does not fit well with any expectations, giving a 
maximum probability for \ 2 of 5% for the best oscillation hypothesis, and the combined probability 
for the shape and number of events is 17%. We conclude that these data favor a neutrino oscillation 
hypothesis, but with unexplained structure in the zenith distribution not easily explained by either 
the statistics or systematics of the experiment. 
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Over the last decade evidence has been growing for 
the possibility of oscillation of atmospheric neutrinos. 
A first anomaly was observed in the ratio of contained 
muon neutrino to electron neutrino interactions in the 
1MB |0| and Kamiokande |^| detectors. In addition, 
the observation of an anomaly in the multi-GeV at- 
mospheric neutrino ratio in Kamiokande suggested spe- 
cific oscillation parameters with large mixing probability 
and Am 2 » 10~ 2 eV 2 ||. Recent results from Super- 
Kamiokande have confirmed the anomaly in the con- 
tained event ratio and also show a strong effect in the 
zenith angle distribution Q] suggesting best fit parame- 
ters of sin 2 29 = 1.0 and Am 2 in the range of a few times 
1CP 3 eV 2 . Also recently, the Soudan 2 detector has con- 
firmed an anomaly in the /i/e ratio of contained events 
using an iron-based detector ^ . Earlier results from the 
Frejus § and NUSEX § detectors are consistent with 
the expected number of contained events though with 
smaller statistics. 

The flux of atmospheric muon neutrinos in the energy 
region from a few GeV up to hundreds of GeV can be 
inferred from measurements of upgoing muons in under- 
ground detectors. If the anomalies in the atmospheric 
neutrinos at lower energy are the result of neutrino oscil- 
lations, then the flux of upgoing muons should be affected 
both in the absolute number of events observed and in 
the shape of the zenith angle distribution, with relatively 
fewer events observed near the zenith than near the hori- 
zontal due to the longer pathlength of neutrinos near the 
zenith. Previous measurements of the upgoing muon flux 
have been made by the Baksan Q , Kamiokande || , 1MB 
|jiof and Frejus || detectors with no claimed discrepancy 
with expectations from calculation. 

The MACRO detector |tl| provides an excellent 
tool for the study of upgoing muons. Its large area 
(76.6 m x 12 m x 9.3 m), fine tracking granularity 
(angular resolution on tracks is between 0.1° and 1.0°), 
good time resolution (about 500 ps), symmetric electron- 
ics with respect to upgoing versus downgoing muons and 
fully-automated analysis permit detailed studies of the 
detector acceptance and possible sources of backgrounds 
to upgoing muons. In addition, the overburden of the 
Gran Sasso Laboratory (minimum rock overburden of 
3150 hg/cm 2 ) is significantly larger than for the locations 
of the previous experiments with the highest statistics on 
upgoing muons (Baksan and 1MB). This provides addi- 
tional shielding against possible sources of background 
induced by down-going muons. 

In our first measurement on upward-going (upgoing) 
muons, we reported on a deficit in the total number of 
observed upgoing muons with respect to expectations and 
also an anomalous zenith angle distribution 

II- Here, 

we report on a data set with much higher statistics |l5| ] 
which retains the same basic features as reported previ- 
ously. In addition, an extensive and exhaustive study has 
been performed on systematic effects in the analysis and 
detector acceptance. 

The upgoing muon data presented here come from 



three running periods and detector configurations: the 
lower half of one supermodule from March 1989 - Novem- 
ber 1991 (1.38 effective live-years), the lower half of 6 
supermodules from December 1992 - June 1993 (0.413 
effective live- years) and the full detector from April 1994 
- December 1997 (2.89 effective live- years). Results from 
the first two periods have already been published jl2| . 

The sign of the direction that muons travel through 
MACRO is determined from the time-of-flight between 
at least two different layers of scintillator counters com- 
bined with the path length of a track reconstructed in 
the streamer tubes. The measured muon velocity is 
calculated with the convention that muons going down 
through the detector will be expected to have 1//3 near 
+1 while muons going up through the detector will be ex- 
pected to have 1/(3 near -1. Several cuts are imposed to 
remove backgrounds caused by radioactivity in near co- 
incidence with muons and showering events which may 
result in bad time reconstruction. The most important 
cut requires that the position of a muon hit in each scin- 
tillator as determined from the timing within the scin- 
tillator counter agrees to within ±70 cm of the position 
indicated by the streamer tube track. 

It has been observed that downgoing muons which pass 
near or through MACRO may produce low-energy, upgo- 
ing particles. These could appear to be neutrino-induced 
upgoing muons if the down-going muon misses the de- 
tector fy}] . This background has been suppressed by 
imposing a cut requiring that each upgoing muon must 
traverse at least 200 g/cm 2 of material in the bottom 
half of the detector. Finally, a large number of nearly 
horizontal (cos 9 > —0.1), but upgoing muons have been 
observed coming from azimuth angles (in local coordi- 
nates) from 30°-50°. This direction corresponds to a cliff 
in the mountain where the overburden is insufficient to 
remove nearly horizontal, downgoing muons which have 
scattered in the mountain and appear as upgoing. We ex- 
clude this region from both our observation and Monte- 
Carlo calculation of the upgoing events. 

Figure H shows the 1//3 distribution for the MACRO 
data from the full detector running (for the older data see 
the equivalent figure in Ref. (l^]). A clear peak of upgo- 
ing muons is evident centered on 1//3 — —1. There are 
398 events in the range —1.25 < 1//3 < —0.75 which we 
define as our upgoing muon sample from this data set. 
We combine these data with the previously published 
data (with 4 additional events due to an updated anal- 
ysis) for a total of 479 upgoing events. Based on events 
outside the upgoing muon peak, we estimate there are 
9 ± 5 background events in the total data set. In ad- 
dition to these events, we estimate that there are 8 ± 3 
events which result from upgoing charged particles pro- 
duced by downgoing muons in the rock near MACRO. 
Finally, it is estimated that 11 ± 4 events are the re- 
sult of interactions of neutrinos in the very bottom layer 
of MACRO scintillator. A statistical subtraction from 
the data is made for these backgrounds prior to calcula- 
tion of the flux. Hence, the observed number of upward, 
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through-going muons integrated over all zenith angles is 
451. 

A Monte Carlo has been used to estimate the expected 
number of upgoing muons. We use the Bartol neutrino 
flux |lj|, which has a systematic uncertainty of ±14%, 
taking into account the agreement with measurements of 
the flux of muons in the atmosphere. We use the Morfin 
and Tung parton set Si jl6) for calculation of the vN 
cross section. These parton distributions were chosen 
based on good agreement of the resulting <jt compared 
to the world average at E v = 100 GeV. We estimate a 
systematic error of ±9% on the upgoing muon flux due to 
uncertainties in a{vN), including low-energy effects 0. 
The energy loss for muons propagating through rock is 
taken from Lohmann et al. |l£j], adjusting the energy loss 
for the average composition of rock in the Gran Sasso. A 
5% systematic uncertainty in the flux of upgoing muons 
results from this calculation due to uncertainties in the 
rock composition and uncertainties of muon energy loss. 
Adding in quadrature all the quoted errors results in a 
total systematic uncertainty of 17% on the expected flux 
which is almost uniform with zenith angle. The expected 
upgoing muon fluxes based on different neutrino fluxes 
p0| pl[ are within 10% of the value presented here. The 
detector has been simulated using GEANT jl9| , and sim- 
ulated events are processed in the same analysis chains 
as the data. An efficiency factor of 0.97 is applied to 
the expected number of events based on various elec- 
tronic efficiencies which have been explicitly measured 
using downgoing muons. 




FIG. 1. Distribution of 1/(5 for all muons in the data 
set taken with the full detector apparatus. A clear peak of 
upgoing muons is evident centered on 1//3 = — 1. The widths 
of the distributions for upgoing and downgoing muons are 
consistent. The shaded part of the distribution is for the 
subset of events where three scintillator layers were hit. 
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FIG. 2. Zenith distribution of flux of upgoing muons with 
energy greater than 1 GeV for data and Monte Carlo for the 
combined MACRO data. The solid curve shows the expecta- 
tion for no oscillations and the shaded region shows the uncer- 
tainty in the expectation. The dashed line shows the predic- 
tion for an oscillated flux with sin 29 = 1 and Am 2 = 0.0025 
eV 2 . 

Care has been taken to ensure a complete simulation 
of the detector acceptance in the Monte Carlo and to 
minimize the systematic uncertainty in the acceptance. 
Comparisons have been made between several different 
analyses and acceptance calculations, including separate 
electronic and data acquisition systems. Studies have 
been made on trigger inefficiencies, background subtrac- 
tion, streamer tube efficiencies, and efficiencies of all data 
quality cuts. Data distributions over many different vari- 
ables (positions of events, azimuth angle, time distribu- 
tions, etc.) have been studied and shown to be consis- 
tent with expectations. The sum (in quadrature) of all 
the systematic errors on the acceptance is ±6% for the 
total number of events. The systematic uncertainty on 
the acceptance for zenith angle bins around the horizon 
is larger than near the vertical due to detector geometry 
effects and smaller statistics for downgoing muons. 

The number of events expected integrated over all 
zenith angles is 612, giving a ratio of the observed num- 
ber of events to the expectation of 0.74 ±0.036(stat) 
±0.046(systematic) ±0.13 (theoretical). The probability 
to obtain a result at least as far from unity as this is 
0.0003 if the Bartol Monte Carlo represents the true 
parent flux of neutrinos. However, taking into account 
the relatively large theoretical uncertainty on the flux 
(mostly on the normalization), the same probability is 
0.14. Hence, there is a low probability that the Bartol 
neutrino flux represents the true flux of upgoing neutri- 
nos at MACRO, but the uncertainty on the normalization 
of this flux makes it difficult to conclude (from this test 
alone) that new physics, such as neutrino oscillations, 
must be responsible for the discrepancy. 
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Figure |2j shows the zenith angle distribution of the 
measured flux of upgoing muons with energy greater than 
1 GeV for all MACRO data compared to the Monte Carlo 
expectation for no oscillations (solid line) and with an 
oscillated flux with sin 2 29 = 1 and Am 2 = 0.0025 eV 2 
(dashed line). The range for the Monte Carlo expecta- 
tion for the unoscillated flux reflects the ±17% systematic 
uncertainty in that prediction. The shape of the angu- 
lar distribution is different than the expectation giving a 
X 2 = 26.1 for 8 degrees of freedom (probability of 0.001 
for a shape at least this different from the expectation) 
for the case of no oscillations but with the number of 
events in the Monte Carlo normalized to the number in 
the data. 




FIG. 3. Probabilities for obtaining the observed MACRO 
results on upgoing muons for — > v T oscillations with 
sin 2 26 — 1.0 and for Am 2 as shown. For the number of 
events, the curve shows the probability of observing a num- 
ber of events which differs from the expectation by at least 
as much as the MACRO data for a given value of Am 2 . For 
the angular distribution, the curve shows the probability to 
observe a distribution which is at least as unlike the expecta- 
tion based on a \ 2 comparison of the shape of the data as a 
function of zenith angle. 

To test oscillation hypotheses, we calculate the inde- 
pendent probability for obtaining the number of events 
observed and the angular distribution for various oscilla- 
tion parameters. Figure || shows the probability of ob- 
taining a number of events which differs from the expec- 
tation by at least as much as the MACRO observation for 
sin 2 29 — 1.0 and various Am 2 for — » v T oscillations. 
(This is a two-sided Gaussian probability.) The expec- 
tation for Am 2 = 0.002 eV 2 agrees with the observed 
number of events. 

The probability of x 2 for obtaining the observed shape 
of the angular distribution has been computed as above 
for oscillation hypotheses and is also presented in Fig. 
|§| for — > v T oscillations. The number of events un- 
der different flux hypotheses is always normalized to the 
observed number of events for this comparison. A max- 
imum probability of 5% is obtained for a distribution at 



least this different from the expectation for — * v T oscil- 
lations. This occurs for sin 2 20 = 1.0 and Am 2 = 0.0025 
eV 2 , but the probability is changed little within a decade 
of Am 2 around this value. However, it is notable that 
the same best value for Am 2 is obtained independently 
from both the angular distribution and the number of 
events. The somewhat low probability for any of these 
hypotheses is the result of the relatively low number of 
events in the region —1.0 < cos# < —0.8 compared to 
the number of events in the region — 0.8 < cos 9 < —0.6. 

Figure ^ shows probability contours for oscillation pa- 
rameters using the combination of probability for the 
number of events and \ 2 of the angular distribution. The 
best- fit point has a probability of 17%. The solid lines 
show the probability contours for 10% and 1% of the best- 
fit value (i.e. 1.7% and 0.17%). The dashed lines show 
the exclusion contours at the 90% and 99% confidence 
levels based on application of the Monte Carlo prescrip- 
tion of reference j25j. The "sensitivity" (not shown) is 
slightly larger than the curve for 10% of P max . The sen- 
sitivity is the 90% contour which would result from the 
preceding prescription if the data and Monte Carlo hap- 
pened to be in perfect agreement at the best-fit point. 
It should be noted that this prescription for producing 
confidence-level intervals assumes that the hypothesis is 
correct. 

Possible systematic effects have been studied and 
shown to be too small to explain the observed anoma- 
lous shape in the zenith distribution. The detector accep- 
tance is best understood (from downgoing muons) near 
the vertical, where the biggest deviation compared to the 
Monte Carlo without oscillations is observed. The data 
from all running periods are consistent in the shape of 
the zenith distribution. We have compared the zenith 
distribution of down-going muons with a Monte Carlo 
expectation based on the known overburden; the two dis- 
tributions agree well. We have compared the measured 
flux of downgoing muons using the same analysis as for 
the upgoing muons and find the result is consistent with 
expectations (see ref. |2(|). The possibility of a water- 
filled cavern below MACRO has been studied, although 
no such caverns are known to exist in the region of the 
Gran Sasso. If all of the region below MACRO were wa- 
ter, a maximum 15% depletion would be observed in the 
flux of upgoing muons. Any realistic water-filled cavern 
would result in a depletion of no more than about 5%. 
For MACRO, we have shown that upgoing charged par- 
ticles produced by downgoing muons contribute a back- 
ground of 2% of the total number of upgoing muons Q . 
This rate could be higher for experiments located in lab- 
oratories with less overburden than the Gran Sasso. 

It has recently been suggested that oscillations between 
and a sterile v could qualitatively produce a shape 
in the zenith distribution of upgoing muons similar to 
that observed by MACRO This would result from a 
matter effect in the center of the Earth. However, due to 
suppressed oscillation amplitude, the current model does 
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not offer a better quantitative agreement with MACRO 
data than the — > v r hypothesis, giving a maximum 
probability of 2%. 

In conclusion, we have reported on a measurement of 
the flux of upgoing muons, produced by neutrinos (with 
< E v >^ 100 GeV) originating in atmospheric cosmic- 
ray showers. The ratio of the number of observed to ex- 
pected events integrated over zenith angles from — 1.0 < 
cosfl < -0.1 is 0.74 ±0.036(stat) ±0.046(systematic) 
±0.13(theoretical). The observed zenith distribution 
does not fit well with the expectation, giving a maximum 
probability for x 2 of only 0.1%. The acceptance of the 
detector has been extensively studied using downgoing 
muons, independent analyses, and Monte-Carlo simula- 
tions. The remaining systematic uncertainties cannot be 
the source of the discrepancies between the data and ex- 
pectations. We have investigated if the anomaly could 
be the result of neutrino oscillations. Both techniques 
independently yield mixing parameters of sin 2 29 = 1.0 
and Am 2 of a few times 10 -3 eV 2 . However, the ob- 
served zenith distribution does not fit well with any ex- 
pectations, giving a maximum probability for \ 2 °f only 
5% for the best oscillation hypothesis. We conclude that 
these data favor a neutrino oscillation hypothesis, but 
with unexplained structure in the zenith distribution not 
easily explained by either the statistics or systematics of 
the experiment. 
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FIG. 4. Probability contours for oscillation parameters 
for — > v T oscillations based on the combined probabili- 
ties of zenith shape and number of events tests. The best-fit 
point has a probability of 17% and iso-probability contours are 
shown for 10% and 1% of this value (i.e. 1.7% and 0.17%). 
The dashed lines show exclusion confidence intervals at the 
90% and 99% levels calculated according to reference |p5f . 
Since the best probability is outside the physical region the 
confidence intervals regions are smaller than the one expected 
from the sensitivity of the experiment. The "sensitivity" con- 
tour (not shown) is slightly larger than that for 10% of P max . 

We are analyzing other topologies of neutrino events. 
We will publish shortly the complementary results from 



semi-contained events with the neutrino interaction 
within the detector (< E v >~ few GeV) and upward 
going stopping muons O], 
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